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We develop a geometric approach to quantify the capability of creating entanglement 
for a general physical interaction acting on two qubits. We use the entanglement measure 
proposed by us for A^-qubit pure states (Phys. Rev. A 77, 062334 (2008)). This geo- 
metric method has the distinct advantage that it gives the experimentally implementable 
criteria to ensure the optimal entanglement production rate without requiring a detailed 
knowledge of the state of the two qubit system. For the production of entanglement in 
practice, we need criteria for optimal entanglement production which can be checked in 
situ without any need to know the state, as experimentally finding out the state of a 
quantum system is generally a formidable task. Further, we use our method to quantify 
the entanglement capacity in higher level and multipartite systems. We quantify the en- 
tanglement capacity for two qutrits and find the maximal entanglement generation rate 
and the corresponding state for the general isotropic interaction between qutrits, using 
the entanglement measure of A^-qudit pure states proposed by us (Phys. Rev. A 80, 
042302 (2009)). Next we quantify the genuine three qubit entanglement capacity for a 
general interaction between qubits. We obtain the maximum entanglement generation 
rate and the corresponding three qubit state for a general isotropic interaction between 
qubits. The state maximizing the entanglement generation rate is of the GHZ class. To 
the best of our knowledge, the entanglement capacities for two qutrit and three qubit 
systems have not been reported earlier. 

PACS numbers: 03.67.Hk, 03.65.Ca, 03.65.Ud 



I. INTRODUCTION 

It is by now well established that entanglement in multipartite quantum systems is a 
physical resource used to perform a variety of information processing tasks [1] as well as 
novel communication protocols [2]. A quantum system evolves to generate entanglement 
provided its parts interact. For such an interaction, the Hamiltonian of the total system 
is not just a sum of the Hamiltonians pertaining to each part (local Hamiltonians). Thus, 
for a bipartite system AB, Hab Ha + Hb but has a term which couples the two 
parts A and B. Together with local operations, this coupling can be used to generate 
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entanglement [3,4,5], to transmit classical and quantum information [4,6,7,8] and more 
generally, to simulate the dynamics of some other Hamiltonian (say H'^^) and thus to 
perform arbitrary unitary gates on the composite space Hab = T^a ® 'Hb [9,10,11]. 

A lot of experimental work is devoted to creating entangled states of quantum systems, 
including those in quantum optics, nuclear magnetic resonance and condensed matter 
physics [12]. Determining the ability of a system to create entangled states provides a 
benchmark of the "quantumness" of the system. Furthermore, such states can ultimately 
be put to some information processing task like superdense coding [13], or quantum tele- 
portation [14]. 

The theory of optimal entanglement generation can be approached in different ways. 
Ref. [3] considers single shot capacities. For two qubit interaction, without any ancilla 
qubits, Ref. [3] presents a closed form expression for the entangling capability and optimal 
protocols by which it can be achieved. In contrast, Ref. [4] considers the asymptotic 
entanglement capacity, allowing the use of ancillary systems and shows that when ancillas 
are allowed, the single shot and asymptotic capacities are in fact the same. However, such 
capacities are difficult to calculate because the ancillary systems may be arbitrarily large. 
In this paper we exclusively deal with the single shot entanglement capacity. Throughout 
this paper, we take h= 1. 

The paper is organized as follows. In section II we deal with entanglement capacity 
for two qubit states while in section III we deal with this problem involving two qutrits. 
In section IV we address the problem of the entanglement capacity involving the genuine 
tripartite entanglement for three qubits. The discussion of the results for two qubit, two 
qutrit and three qubit cases is included separately in sections II, III and IV respectively. 

II. THE TWO QUBIT CASE 

We develop a geometric approach to calculate the entanglement capacity of any two 
qubit system interacting via a Hamiltonian which is locally equivalent to 

Hi = i_iia^ ® o-f + H2(T2 ® o-f + 1^3(^3 ® (rf . (1) 

Here > ^2 ^ A^a- In this section we omit qubit and system identifiers A, B and AB. 
We define the single shot entanglement capacity by 

F"*"^ = max lim^ ^- . (2) 

|V>G«7t-+0 t 

The Hamiltonian Hi in Eq. (2) is given by Eq. (1). E{\ip)) in Eq. (2) stands for the two 
qubit pure state entanglement measure given by us and is shown to have all the essential 
(as well as many desirable, e.g., superadditivity and continuity) properties expected of a 
good entanglement measure [15]. For a A"-qubit pure state 

^(l^)) = lir(^^ll-i 
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where | IT*-^-* 1 1 is the Hilbert-Schmidt (Euchdean) norm of the N way array X^-^) occurring 
in the Bloch representation of [15,16]. 

The scenario we address, is as follows [3]. The idea is to supplement the interaction 
Hamiltonian Hj with appropriate local unitary operations in such a way that the state 
of the qubits at any time t is precisely \iJjE{t)), for which the increase of entanglement is 
optimal. In order to construct such a procedure, we consider the evolution given by Hj 
to proceed in very small time steps St. Let us also assume that the qubits are initially 
disentangled. Using local operations, we can always prepare the state \ipo) that is, the 
product state which most efficiently becomes entangled under the action of Hj. After a 
time step dt, the state will change and its entanglement will increase to 6E. Then, we 
use (fast) local unitary operations to transform the new state of the qubits into the state 
IV'^e) for which F is optimal. Note that this is always possible, since for qubits all states 
with the same value of E, say 5E, are connected by local unitary transformations. By 
proceeding in the same way after every time step, and taking the continuous time limit 
St — 7> 0, we obtain that the state of the qubits at time t is always the optimal one, \ipE{t))- 
Obviously, in an experimental realization, this procedure requires that we can apply the 
appropriate local transformations in times which are short compared to the typical time 
scale thj associated with Hj, th = {cmax — emin)~^ , where Cmax and Cmm are the maximum 
and minimum eigenvalues of Hj. Note that Eq. (2) defines the entanglement capacity as 
the maximum achievable entanglement rate for a given system with given interactions. 
We are also interested in finding the state \ipmax) for which the entanglement rate is 
maximum, (denoted by F"^"^ in Eq. (2)). 

We consider two qubits interacting via the Hamiltonian Hi in Eq. (1), which represents 
general interaction between two qubits [3]. First we find the entanglement rate F given 

by 



lim 



Here is given by a general two qubit state in the Bloch representation [15,16], 

P= |^)(^| = ^ I /«)/ + 5]]rfcafc®/ + ^s,/(8)(7, + ^Tfe,(7fe(8)(7H , (4) 

\ k I k,l / 

where a^^i, /c, / = 1, 2, 3 are the Pauli operators. We denote by T = [r^j] the correlation 
matrix occurring in the last term of Eq. (4) . Tij are defined by 

Tij = Tr{ai (g) ajp) = {ip\ai ® crj#)- (5) 

Tfe and si are the components of the Bloch vectors [16] of the reduced density operators 
Pa and ps respectively, given by 

rk = Tr{akPA) = i^Wk'^Il'ip), (6a) 



3 



Si = Tr{aipB) = 
We define the entanglement of tlie state l^^) as [15] 

^(IV')) = lir||-i, 



(66) 



(7) 



wliere ||T|| = \jY^ij=i "^ij Euclidean norm of T. For two qubits, this measure is 

related to concurrence [15] and hence to the Von Neumann entropy of the reduced density 
matrix. 

After finding T we maximize it, using a simple geometric argument. It is heartening 
to see that the scenario described above emerges naturally out of this geometric method. 
The entanglement rate F is given by (See Eq. (3)) 

dE_d\\r\\_ 1 ^ 



dt 



dt 



\\r\\ 



with Tij given by Eq. (5). We evaluate fy as follows. 



dr, 



"Tij — 



d 



dp^ 



We now use the equation of motion , 

where the Hamiltonian Hj is defined via Eq. (1), to get [17], 



dr, 



dt 



-iTr{(Ti ® aj[Hi, p\) = iTr{Hi[ai (g) aj, p]). 



Substituting p from Eq. (4) and using the commutation relations [17] 

[ai (g) aj, ak ® ai] = ^[ai, ak] {aj, ct/} + ^{cTi, ak} ® [aj, ai], 

we get, using [uj, aj] = 2iSijk(Jk, {uj, aj} = 2Sij and the expression of Hj in Eq. (1), after 
a bit of algebra. 



dr. 



dt 



= -2 



l,n 



k.n 



This gives 



fij — 2 



j,k,n 



i.k.n 



Thus we get, for the entanglement rate F, 
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] /^n- (8) 

n 

Here -r„ and Tn-. are, respectively, the n th column and row vectors of the correlation 
matrix T = [rij] . 

The entanglement generation rate T expressed in Eq. (8) is obtained via the temporal 
evolution of the initial state by the interaction Hamiltonian Hj. This expression for F does 
not depend on any local unitary transformation applied to a qubit. Following the general 
scenario described above, (see the third paragraph of this section), we now lock on to an 
instant of time and apply the local unitary transformations to qubits, in order to find the 
conditions for optimal F and the corresponding two qubit state \iPe)- The experimental 
meaning of this sentence is described as a part of the scenario above. In the geometrical 
approach we have adopted, local unitary transformations amount to rotations of vectors 
in Eq. (8), which are the vectors in the Bloch space of individual qubits. Wc expect the 
entanglement to remain unultered by the local unitaries, which turns out to be the case. 
The entanglement measure in Eq.(7) is not affected by local unitaries, as proved in [15]. 

Obviously, F will be maximum if the components of the vector products occurring in 
Eq. (8) are replaced by the magnitudes of these vector products and the factors in these 
products are mutually perpendicular. Geometrically, this means that the vector products 
themselves are in the directions of the components occurring in Eq. (8) with the other 
two orthogonal components zero. For example, {fx t-i) is along its first component, i.e. 
along X axis, with its y and z components zero. Thus, in order to maximize the first term 
in Eq. (8), namely, 

^{rX r:n)nl^n = X -^1)1^1 + (f X -^2)2/^2 + (r X ^3)3/13, 

n 

we must have vectors {fx f-i), {fx 7^2) and {fx 7^3) along x,y, z axes respectively. This 
can be done only when one of the vector products is zero. Since /^i > /i2 > A^s, we choose 
(fx 7^3) = 0. Given the vector (fx r.i) along the x axis and the vector (fx 7^2) along the 
y axis, we can choose f to be along the z axis and vectors -ti and 7^2 along the y and x 
axes respectively. In exactly the same way, maximization of the second term in Eq. (8), 
^^(s X r„:)„yLi„, makes the vector s along the z axis and vectors ti- and f2; along the y 
and X axes respectively. Writing explicitly the components of the vector products in the 
expression for F (Eq. (8)) and putting ri_2 = = Si^2 we get, 

o 

((-r3T2l - 53X12)//! + (r3Ti2 + S3T2l)At2)- 



Since we are dealing with the two qubit pure states we have ||f|| = [18], so that 
r3 = ±53. Choosing 

r3 = -S3 (9) 
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we get, 

2 

r = ||:^r3(Ti2 - r2i)(Aii + /X2). 

The expression (ti2 — T21) becomes maximum when 

ri2 = -r2i- (10) 

Finally, we note that this maximization procedure docs not change ||T|| = v^^^jF^nlF 
and hence the entanglement value given by Eq.(7). Further, choosing the cross products 
along their components appearing in Eq. (8) corresponds to the rotations in Bloch space, 
generating local unitaries on the system. Therefore, the maximum of F over the states 
with same entanglement, that is, F^, is given by 

4 

Fe = TT7piTr3ri2(Ati + (J't)- (11) 



To get the state \iPe) corresponding to F^, we seek the state satisfying conditions Eq. 
(9) and Eq. (10). We start with the general state \iJj) = j=oCij\ij) and calculate T12 
and T21. In order to satisfy Eq.(lO), the state IV') should be 

IVe) = |coi||01)+i|cio||10) ; |coi|'+|cio|' = l, (12) 

which is the same as \iPe) obtained in Ref [3] if we identify |coi| = ^/p. Further, we can 
write Te (Eq. (11)) as the product of two factors 

= f{p)hmax 

with 
and 

m - (13) 



To get /(p) as a function of p, we calculate r^, and ||T|| using the state \iPe) (Eq- 
(12)) so that 

f( ) = 113^2 ^ 8(l-2p)v/p(l-p) 

lirii vi + 8p(i-p) 

Fig.(la) depicts this f{p) verses p, while Fig. (lb) plots the analogous f{p) obtained us- 
ing Von Neumann entropy as the entanglement measure. Note that f{p) and hence F^ 
vanishes for the maximally entangled state (p — \) about which it is antisymmetric 
/(| + x) = — /(| — x). We see that, as p increases from to |, F^; > makes the en- 
tanglement increase, until is maximal at p = |, after which F < 0, making entanglement 
decrease to zero as p approaches 1. 
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Fig. la 




Fig. lb 

Fig.l: (a) f{p) vs. p for entanglement measure in Eq. (7) and (b) fvN{p) vs. p for Von 
Neumann entropy of the reduced density matrix (see text) . 
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Thus we see that, F^; is the product of the function which depends only on the state, 
(via p) and the factor hmax which depends only on the interaction strengths fii and /i2, that 
is, on the interaction Hamiltonian. Note that hmax is independent of the entanglement 
measure. The form of f{p) for the entanglement measure in Eq.(7) and that for the Von 
Neumann entropy, (Fig. (lb)) also turns out to be the same. To get F"^"^ we have to find 
Po at which f{p) is maximum. To do this, we invoke the relation between fvN{p), which is 
the analog of f{p) in Eq. (13) obtained via Von Neumann entropy of the reduced density 
operator of the two qubit pure state [3] and f{p) in Eq. (13) obtained via E{\iIj)) in Eq. 
(7). This is 



where Evn is the Von Neumann entropy of the reduced density operator and E is given 
by Eq. (7). Maximizing the RHS of Eq. (JH]) we get po « 0.0832217 and F""''^ ^ 1.9123. 
The state \ipmax) corresponding to T"^°-^ is the state \iIje) with p = pq. 

From the definitions of and rfc(Eq. (5) and (6a)), these quantities are averages of 
the Pauli operators in a state {ip), which can be obtained using experimentally measured 
values of the corresponding operators on two qubits. Therefore, the geometric method 
presented here has the advantage that the conditions for Te, Eqs. (9) and (10), can he 
tested experimentally giving us an experimental way to check out whether the system has 
reached the state \iPe)- The value of F^; can also be experimentally estimated via Eq. 
(11), for given /ii + /i2- Further, the function f{p) can be estimated experimentally via 
Eq. (13) as the system evolves, under the given Hamiltonian, toward \ipmax)i or under the 
local unitaries toward \iPe)- These experimental estimations can he carried out without a 
detailed a priori knowledge of the quantum state at any time during the evolution of the 
two quhit system. These facts can be of great advantage in a practical implementation of 
any scheme to entangle two qubits interacting via some Hamiltonian or quantum gates 
[19]. For the production of entanglement in practice, we need criteria for optimal entan- 
glement production which can be checked in situ without any need to know the state, 
as experimentally finding out the state of a quantum system is generally a formidable 
task. We note that, in order to achieve such an experimental determination of optimal 
entanglement production rate using the model in ref [3], we have to experimentally obtain 
the values of the Schmidt coefficients of the evolving two qubit state, which requires the 
experimental determination of the two qubit state itself. This requires more experimen- 
tal effort and resources {d + 1 different joint measurements, d = dimension of the joint 
Hilbert space [20]) as compared to measuring the quantities in Eqs. (9) and (10), which 
are simply the average values of the Pauli operators in the state. 



The entanglement measure in Eq. (7) can be generalized to the qudit pure states 
which satisfies all the essential (and many desirable, e.g., superadditivity and continuity) 




(14) 



III. THE TWO QUTRIT CASE 
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properties expected of a good entanglement measure [21]. Therefore, we can use it to 
obtain the entanglement generation rates for the multipartite d level systems. Here we 
find the entanglement generation rate for two qutrits (labeled A and B) interacting via 
the Hamiltonian 

8 

p=l 

where the interaction strengths satisfying /i^ > for k < l,k,l = 1, . . . , 8. Here 

Ap, p = 1, . . . , 8 are the generators of the SU{3) group satisfying Tr{XiXj) = 26ij and are 
characterized by the structure constants of the corresponding Lie algebra, fijk and gijk, 
which are, respectively, completely antisymmetric and completely symmetric. 



Other useful relations are 



iifjki = Tr{[\i, \k]\i) 
4:gikp = Tr{{Xi, Xk}Xp). 



(16) 

(17) 
(18) 



We give here the generators of SU{3) in the |1), |2), |3) basis [17] to be used below. 
Ai = |l)(2| + |2)(l| 
A2 = -.(|l)(2|-|2)(l|) 
A3 = |l)(l|-|2)(2| 
A4=|l)(3| + |3)(l| 
A5 = -^(|l)(3|-|3)(l|) 
A6 = |2)(3| + |3)(2| 
A, = -.(|2)(3|-|3)(2|) 
A8 = ^(|l)(l| + |2)(2|-2|3)(3|). 

The action of these generators on the basis states { 1 1) , |2) , |3) } is given by the following. 

Ai|l) = |2), Ai|2) = |1), Ai|3) =0 
A2II) = z|2), A2I2) = -z|l), A2I3) = 
A3|l) = |l), A3|2) = -|2), A3|3) = 
A4|l) = |3), A4|2) = 0, A4|3) = |l) 
A5|l)=^|3), Ai|2) = 0, Ai|3) = -z|l) 
A6|l) = 0, A6|2) = |3), A6|3) = |2) 
A7|l) = 0, Ar|2)=z|3), Xj\3) = -t\2) 
A8|l) = ^|l), A8|2) = ^|2), A8|3) = -^|3). 

We use these equations to get the vectors A"^ and A"^ in whose components are the 
averages Af = (?/'|Aj (g) = Tr{Xi Ip), i = 1,...,8 and Af = {ipll Xilijj) = 
Tr{I Xip), i = 1, . . . ,8 respectively, where {ip), (p = \ip){ip\) is a two qutrit pure state 
(see Eq. fl20|) and the discussion following it). 



The pure state entanglement for two qutrits is given by [21], 

^(l^)) = lir||-3 



where ||T|| = \jYliij=i'^ij the Euchdean norm of T. The general two qutrit pure state 
p has the following Bloch representation. 



P 
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Here {\^'^) = Tr(pA,B^^'^) with A"^'^ and pa,b (the reduced density operator) apply to 
the qutrit A and B respectively, while r^i = (9/4)Tr(A^ ® Af p). The definitions of ffc; and 
r are 

1 



lirii 



7 ^ "^ij^iji 



fij = iTr{H[Xi (g) Aj,p]), 

where we have used the Heisenberg equation of motion as in the two qubit case. Using 
Eq. ( IT6l) . ( 1T7|) . ( IT8l) and the elements of the tensors fijk and gijk in [17], we get, after 
some algebra, the following expression for F 

r = -3 (-^) Yl ^'pf^^p (""^p^^ + ^p^^n ■ (19) 

^" k,p,l=l 

Expanding the sum in Eq. (|T9|) and rearranging, we get. 



r = -3 (j^^ ^a(5) 5^/i,[(^, X A^), + {t,., X A^),], (20) 
^" s pes 

where S runs over the triplets 

(1, 4, 7), (2, 1, 6), (3, 1,5), (3, 2, 4), (2, 5, 7), (3, 7, 6), (5, 4, 6), (3, 6, 8), (2, 5, 8) 

anda(^) has values 1, 1/2, 1/2, 1/2, 1/2, 1/2, 1/2, v^/2, ^3/2 respectively for these triplets. 
Tp and Tp: are the vectors in with p G 5* where 5* is one of the above triplets and the 
index : varies over a given S for fixed p. X"^'^ are vectors in respectively comprising 
the components of A"^'^ indexed by one of the triplets 5*. There are in all 54 terms in 
Eq. ( l20i) . Unfortunately, all these terms are coupled and a simple geometrical procedure 
to maximize F, as in the two qubit case, seems very difficult. However, it is straight- 
forward to maximize F numerically over the coefficients Cij, i,j = 0, 1,2, by expressing 
all the terms in the expression for F (Eq. (l20ll ) as averages in the general two qutrit 
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state = i,j = 0,1,2. We can carry out the numerical maximization for 

the general Hamiltonian in Eq. ( !T5|l . where the strengths of interaction /i^ have different 
values. In that case, T does not have the simple structure analogous to F = f{p)hmax as 
in the two qubit case. Therefore, we assume isotropic interactions so that all interaction 
strengths are equal to a common value /z. In this case, F has a simple form 

F hjYiaxifiiji iyj 0, 1,2)^. 

Therefore, we maximize F assuming the isotropic interactions. The result is 

ymax _ 3 99495^ 

and the corresponding (normalized) state is given by 

Coo = -0.28317 + i0.148948; coi = -0.433055 + iO.382479; C02 = -0.117778 + iO.274948 

cio = 0.0625717 - i0.144584; cn = 0.102783 - i0.0787094; cu = -0.340939 - z0.324717 

C20 = 0.25066 - i0.167261; C21 = 0.0344755 - iO.244282; C22 = 0.227159 - i0.088347. 

After converting this state to the Schmidt canonical form we get the state giving the max- 
imal entanglement generating rate for two qutrits, under a general isotropic interaction, 
as 

li^max) = 0.884297|00) + 0.448838|11) + 0.128697|22). 

We find that Edipmax)) = 0.677882. This shows that, in order to increase the entan- 
glement of a two qutrit system in an optimal way, it is better to start with an initially 
entangled state rather than a product state, at least when all the interaction strengths in 
Hi (Eq. (fT5l) ) are equal. We also note that the optimal entanglement E{\ipmax)) is inde- 
pendent of Hj, provided, again, that all interaction strengths in Hj (Eq. (fT5l) ) are equal. 



IV. THE THREE QUBIT CASE 

We now deal with the problem of entanglement generation capacity for three qubits. 
We emphasize that this is the capacity to generate genuine three qubit entanglement and 
not the bipartite entanglement between any two parts of the three qubit system. In this 
case also, the entanglement measure given by Eq. (7) can be used as this entanglement 
measure applies to A^-qubit pure states and has all the essential (and many desirable, eg 
superadditivity and continuity) properties expected of a good entanglement measure [15]. 
For the three qubit case, r in Eq. (7) is the three qubit correlation tensor appearing in 
the Bloch representation of the state. Here r is a three way array while for two qubits r 
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was a matrix. The Bloch representation of a general three qubit pure state is 
p = |^)(^/,| = ^{l ® I ® I + ^riai® I ® I + ^SnI ® (Tn® I + ^qml ® I ® (7m + 

I n m 

In Im nm 

+ ^ rimn<yi (Tn (g) (T^) . (21) 

Imn 

Here r — [rijk] is a three way array while t^"^ — [t\j^] are matrices. The definitions of 
various symbols in p are as follows. 

n = Tr{atpA)^Tr{at®I®Ip) 
Sr, = Tr{a^pB)^Tr{I®a^®I p) 
Qm = Triagpc)^Tr{l0l<^agp), 



C = TT{at®a^pAB) = TT{at®a^®I^p) 
tf^ = Tr{at®agpAc)^Tr{at®I^®agp) 
tnm = rr(a^®a^pBc)=rr(7^®<j^®<7^p), 

where Pa,b,c and Pab,ac,bc ai'e the appropriate reduced density operators. We consider 
the general interaction between qubits which can be reduced by the singular value decom- 
position to the Hamiltonian 

Hi — Hab + Hac + Hbc (22) 

where 



tC 



Hab = E^^V®^s^^^' 
Hbc = i^^^f^I' 

s=l 

Hac = (23) 



It is helpful to imagine that the three spins are at the vortices of a triangle. If they are 
arranged on a line, we expect on physical grounds that one of the terms can be neglected 
in comparison with the other two, as it gives the next nearest neighbor interaction. For 
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all subsystems we have /ii > /i2 > /^s- Using the definition of the entanglement generation 
rate F in Eq. (3) and the definition of the entanglement measure in Eq. (7) we get, 



1 ^ 

^ ^ TItTI ^ 'Tijk'T'ijk- 



i,j,k=l 



Using the Heisenberg equation of motion, 



where the Hamiltonian Hj is defined via Eq.s (I22l) . (l23|) . we get 

fijk = iTr{Hi[ai ® cij a^, p]). 
We now use the commutator identity 

[A®B®C,D®E®F] = h^[A,D]®[B,E]® [C, F] + [A, D] ® {B, E} ® {C, F} 

+{A, D} ® [B, E] ® {C, F} + {A, D} ® {5, E} ® [C, F]) 
and the definitions of p and Hi in Eq.s (!2T|) and (!22l) respectively to get, 

3 3 



ij/j e^ij'fc + Pk 2^ ^ifc' ^jfc'fc + 

j'=l k'=l 

3 3 

2^i^ii' ^ki'j + fJ-i ; ^I'jv ^ki'i + 



/'=1 

3 



/'=1 
3 



"•"Pi ^ ^k'k ^jk'i -r Pj Ij'fc 

j'=l 



(24) 



k'=l 



where e s are the Levi-Civita symbols. Substitution of Eq. ( 124|) in the expression for F 
gives, 



F 



-2 

M 



k,s=l 



i,s=l 
3 



AC 



(25) 
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where r-sk = [Tisk,T2sk,T3sk]'^ , for example, is a vector in for fixed s and k. Similarly, 
t'.'i^ and tj. are the kth. column and the jth row vectors of the matrix t^"\ The expression 
for the entanglement generation rate F for three qubits (Eq. ( 125|) ) has 54 coupled terms, 
each term being a component of the cross product of two vectors. A geometric argument 
to maximize F, as in the two qubit case, seems to be very difficult. However, it is quite 
straightforward to maximize F numerically, by writing the elements of the three way array 
T and the matrices t^' -* as the appropriate averages in the general three qubit state 

7 

where i labeling the product basis ket \i) is the binary representation of the index i. We 
can numerically optimize F for the general Hamiltonian in Eq. fl2^ . However, for the 
general case, where the interaction is anisotropic, that is, the strengths of interaction /i^"'' 
have different values, F does not have the simple structure F = f{p)hmax as in the two 
qubit case. Therefore, we assume isotropic interactions so that all interaction strengths 
are equal to a common value fi. In this case, after evaluating all the terms in Eq. f l25|) in 
the state given by Eq. f l2B]) . F can be written as 

F = /l(co, . . . ,C7)/i. 

After the numerical optimization of F as a function of Cj, i = 0,1, . . . , 7, we get. 

The (normalized) state corresponding to this T"^"-^ is given by 

I'ipmax) = (0.033768 - i0.168758|000)) + (0.574022 - i0.0709471)|001) + 

+ (0.0218412 - i0.111565)|010) + (0.672021 - i0.0754116)|011) + 

+ (-0.0603488 + i0.172566)|100) + (-0.0051137 - i0.183831)|101) + 

+ (0.0556843 + i0.151888)|110) + (0.0700719 - 20.259423)|111). 

This state has the following Acin canonical form, expressed by the two fold degenerate 
sets of entanglement parameters [22]. 

\^+) = 0.610291|000) + 0.67402 exp(z2.51395) 1 100) + 0.394893|101) + 

+0.110357|110) + 0.0715772|111), 

or, 

\^p^) = 0.329873|000) + 0.546087 exp(z0.402558)| 100) + 0.730583|101) + 
+0.20417|110) + 0.132424|111). 
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We see that the state with the maximal entanglement generation rate F"^"^^ belongs to 
the GHZ class. Further, we find that Ediprnax)) = 0.258918. This means that, given the 
isotropic interaction, it is beneficial to start with an entangled three qubit state for optimal 
entanglement generation. Also, we note that the optimal entanglement is independent of 
Hj, provided the corresponding interaction is isotropic. 

Thus we see that, for three qubits, the geometric method based on the entanglement 
measure given by Eq. (7) can be numerically implemented to get the state with maximal 
entanglement generation rate. This program can be carried out for the general interac- 
tion Hamiltonian in Eq. (122|) although we have restricted to the isotropic interactions. 
This procedure can be suitably carried out in a laboratory using quantum circuits. Every 
quantum circuit acts unitarily on a quantum state and we can always find a Hamiltonian 
corresponding to such a circuit [9, 23]. On the other hand, given a (interaction) Hamilto- 
nian for a three qubit system, we may construct a circuit implementing the corresponding 
evolution using universal quantum gates. We note that, for the isotropic interaction, the 
maximal entanglement generation rate F"*"^ is proportional to the interaction strength 
H and the corresponding state is independent of as in the two qubit case. When the 
interactions are anisotropic, the scenario for two qubits does not apply to the three qubit 
case, as r does not factor into the product of a state dependent function and an expression 
involving only the interaction strengths. Even when the interactions are isotropic, we do 
not know the explicit form of such a state dependent function. In other words, we do not 
know whether it is possible to separately account for the contribution due to the state 
and that due to the interactions. Thus a general procedure for the maximization of the 
entanglement generation rate for the higher dimensional and multipartite systems still 
seems to be an open question. These observations ensue from the fact that the terms in 
the expression for F could not be decoupled. This difficulty seems to be generic, as it may 
be a consequence of the difficulties in the geometric interpretation of the Bloch space for 
the multipartite and higher dimensional systems [24]. All the remarks in this paragraph 
apply to the two qutrit case as well. 
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